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Abstract: We study brane inflation in flux compactifications at the nonlinear level, 
solving the D3-brane DBI-equations of motion including backreaction from the D3- 
D7 brane potential and from perturbations of the warp factor. We first numerically 
compute the exact functional form of the Kahler modulus valid on the entire su- 
pergravity background and obtain a two-field potential along radial and harmonic 
directions. We find that a valid perturbative expansion on the entire supergravity 
background with the Kahler modulus integrated out adiabatically in DBI inflation 
requires hierarchies of scales that determine compactification parameters different 
from those typical in slow-roll models. Our numerical results then show that the 
DBI inflationary solutions are quite robust against these nonlinear corrections. 
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1. Introduction 

The isotropy and homogeneity of the universe suggests that in very early times after 
the big bang the universe must have gone through a period of rapid exponential ex- 
pansion called infiation. The quest for deriving infiation from a fundamental theory 
has led to the construction of infiationary models within the framework of string 
theory. For embedding infiation into string theory, one has to specify a string com- 
pactification whose low-energy effective theory contains a suitable inflaton field and 
potential. The many moduli fields that arise from Calabi-Yau compactifications can 
provide a large number of candidate scalar fields in the four-dimensional theory and 
any of these scalar fields may play the role of the inflaton field. However, moduli 
fields are generically either massless, or have a potential with runaway behaviour, 
which makes their interpretation as inflatons rather difficult. In addition to this, the 
computation of the effective potential in terms of all these scalar fields and degrees 
of freedom is a highly nontrivial task. Nevertheless, by considering the systematics 
of fiux compactifications [1] and nonperturbative effects [2] in string theory, it is pos- 
sible to reduce the degrees of freedom and stabilize all the moduli fields as required 
for a suitable model of infiation. 

One particular way of constructing string infiationary models has been within 
the 'KKLT' framework, [2], according to which the older idea of brane infiation [3-6] 
is realised via the motion of a brane in the compact extra dimensions [7]. It was 
shown in [7] that the potential of the D3-brane in a Calabi-Yau fiux compactification 
of type IIB theory [1] with the Kahler moduli stabilized according to [2] receives 
corrections that make the infiaton potential too steep for slow-roll infiation. It was 



- 1 - 



also conjectured in [7] that further compactification corrections to the inflaton poten- 
tial may lead to a finely tuned cancellation of correction terms, so that the infiaton 
potential can be made flat. Major progress in this direction was made when these 
corrections were computed from ultraviolet deformations of the warped throat ge- 
ometry and applied to explicit brane inflation models [8-22] (see also the reviews 
[23-29]). In these models the flattening of the inflaton potential is analysed under 
flne tuning and vast scanning of the compactiflcation parameters. In particular, in 
the 'delicate universe' models studied in [12, 14, 15] it was shown that the inflaton 
potential can be made flat only in small region around an inflection point; elsewhere 
the potential remains steep. However, whether or not the potential can be made flat 
to realise slow-roll conditions, inflation can occur by the DBI effect with high speed 
and steep potentials [30, 31]. But the usual approach taken by the DBI brane infla- 
tion models in the literature is to consider brane motion in an imaginary self-dual 
(ISD) compactiflcation in which the effects of moduli stabilization and backreaction 
are not included. 

Any successful brane inflation model has to include the effects of compactiflcation 
and moduli stabilization, which induce important corrections to the inflaton action. 
One particular correction from moduli stabilization is the departure of the ISD solu- 
tion which induces harmonic dependence in the action of the probe D3-brane, [13-18], 
which has been studied to some extend in [19-21] and analysed in more detail in [32]. 
In previous work [32], we solved the DBI brane equations of motion in the warped 
deformed conifold [33] with harmonic dependent corrections from linearized pertur- 
bations around the ISD solution. We showed that just as angular motion increases 
the number of e- foldings (spinflation) [34-37] , having additional angular dependence 
from linearized corrections also increases the number of e-foldings. However, this 
line of analysis considered multifleld effects in D-brane inflation from only linearized 
perturbations around the ISD supergravity solution. In this paper, we extend our 
previous analysis [32] by including further harmonic dependent corrections from non- 
linear perturbations around the ISD solution, which also contribute to the inflaton 
action. In the noncompact limit, non-linear perturbations are dominated by imagi- 
nary anti-self-dual (lASD) fluxes sourced by moduli stabilizing wrapped D7-branes 
and the flux induced potential for the probe D3-brane in ten-dimensional super- 
gravity equals the nonperturbatively-generated D3-D7 potential in four-dimensional 
supergravity [18]. 

One particular motivation of extending our previous analysis [32] by including 
such nonperturbative corrections is the possibility of increasing the amount of infla- 
tion and decreasing the level of non-Gaussianity due to backreaction effects [17]. The 
backreaction on the mobile D3-brane is sourced by itself [13]. The D3-brane in a flux 
compactiflcation containing holomorphically embedded D7-branes corrects the warp 
factor which in turn corrects the warped volume of four-cycles wrapped by D7-branes. 
This then corrects the D3-D7 potential causing the backreaction on the D3-brane. 
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The perturbations of the warp factor are given by the Green's function and correct 
the 7DBi-factor which controls the level of non-Gaussianity. The Green's function 
solves the noncompact supergravity equation of motion describing the perturbations 
around the ISD solution sourced by lASD fluxes and may be expanded in an infinite 
set of eigenstates of the Laplacian containing harmonic dependent hypergeometric 
functions [38, 39]. 

In this paper we solve the D3-brane equations of motion from the DBI action 
in the warped deformed conifold [33] with harmonic dependence from the leading 
correction terms of the warp factor and the D3-D7 potential. To solve the brane 
equations of motion, we note that the D3-brane potential including nonperturbative 
corrections depends on the functional form of the Kahler modulus and that of the 
D7-brane embedding. For the Kuperstein embedding of D7-branes [40] and a simple 
choice of an harmonic dependent trajectory on the deformed conifold, we integrate 
out the Kahler moduhis and reduce the multifield D3-brane potential to a simple 
two-field potential depending on one radial and one harmonic direction of the coni- 
cal geometry. In our numerical integrations we integrate out the Kahler modulus by 
computing its exact functional form valid on both the infrared (IR) and ultraviolet 
(UV) regions of the supergravity background. We find that computing the Kahler 
modulus within the adiabatic approximation in DBI infiation requires certain hierar- 
chies of scales which determine the set of compactification parameters different from 
those in slow-roll models. Integrating the brane equations of motion for the consistent 
choice of parameters with the numerically computed Kahler modulus shows that the 
DBI inflationary solutions are quite robust against non-linear harmonic dependent 
corrections from perturbations of the warp factor and the D3-D7 brane potential. 

The paper is organized as follows. In section 2, we outline our supergravity 
background and specialize to the warped deformed conifold. In section 3, we consider 
the general set up for brane inflation in the warped deformed conifold deriving the 
explicit form of our D3-brane potential and its relating DBI brane equations of 
motion. In section 4, we present our numerical results. We first determine the 
relevant parameter space and compute the Kahler modulus. We then solve the DBI 
brane equations of motion. In section 5, we conclude with a brief discussion. 



2. Type IIB supergravity 

The setting that we choose for our brane infiationary analysis is the Calabi-Yau fiux 
compactification of type IIB theory containing a warped throat region generated by 
fiuxes [1]. We therefore consider backgrounds in low-energy IIB supergravity, which 
in the Einstein frame can be represented by the action 
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Here 5*100 stands for localized contributions from D-brane and orientifold planes; 
T = Co + ie~^'^ is the axion-dilaton field; G3 = F3 — tH-^ is the combination of 
the R-R and NS-NS three-form fluxes F3 = dC2 and = dB2 (with C2 and 
B2 being the R-R and NS-NN two- forms, respectively); C4 is the R-R four- form; 
F5 — dC^ — ^€2 A + \B2 a F3 is the self-dual five-form, F5 = ^10-^5, with 
denoting the ten-dimensional Hodge-star operator; T^io is the ten-dimensional Ricci- 
scalar; k\q — \{2T^y a'* is the ten-dimensional gravitational coupling. 



In a flux compactification we may take the line element of the form 



dsl,^h-'/\y)dsl-h'/\y)dsl 



(2.2) 



where h is the warp factor depending only on the internal coordinates y := y"", dsl 
is the metric on the internal manifold and ds^ is the four-dimensional flat metric, 
which will be the FRW metric for cosmological analysis. 

Following [1] , we take the self-dual five- form to be given by 



[1 + * 
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da{y) A dx^ A dx^ A dx^ A dx^ 



(2.3) 



in the Poincare invariant case, where a{y) is a function of the internal coordinates. 
The Einstein equations and five-form Bianchi identity then imply 



A(o)$± = ^1 Gil' + /i|V$±| + 7^4 + local 



(2.4) 



where A(o) is the Laplacian with respect to the six-dimensional unperturbed Calabi- 
Yau metric gmn, and 

The equation of motion for the three-form flux is 



— ± a. 
h 



(2.5) 



where by definition 



i dr , . - , 
dA + A (A + A) = 0, 
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For G- = 0, i.e. ^6^3 = iGs , the flux G3 is imaginary self-dual (ISD), the back- 
ground metric is Calabi-Yau with the five-form flux given by a = h^^ [1]. The 
background that satisfies these conditions is called ISD. The specific example of such 
a background that we are interested in is the warped deformed conifold [33] . The de- 
formed conifold is a noncompact and nonsingular Calabi-Yau three-fold in defined 
by the following constraint equation: 

Ei'Ar = e'. (2.8) 

A=l 

where {z^, A = 1,2, 3, 4} represent the local complex coordinates in C^, and e is the 
deformation parameter which can be made real by phase rotation. For vanishing 
e, Eq. ( p. 81 ) gives the singular conifold and describes a cone over a five-dimensional 
Einstein manifold X5. For us the nonsingular limit is relevant in which X5 is the 
[SU{2) X S'f/(2)]/ U{1) coset space T^'^ of topology S"^ x parametrised by a set 
of five Euler angles = {6i,(pi,ilj} with < 6i < it, < < 27r, < ip < An 
{i = 1,2), and the would-be singularity at the tip, r = 0, is replaced by a blown-up 
5''^ of T^'^ amounting to the deformation measured by e. The base of the cone can 
be parametrised by the coordinates yi in a standard way [41] 



y^ = — { cos — cos ^e2^^i+^2+^) _ _i, ^ -^[^,+^2-^) /2.9) 

= ' /^cos-cos-et(^i+^2+'^) + sin-sin-e-^(^^+^^-^)Y (2.10) 
v2 \ 2 2 2 2 J 

y3 = -^( cos ^ sin ^e5(^-^^+'^) + sin ^ cos ^-le'^(^^-^'+A , (2.11) 
v2 \ 2 2 2 2 J 

= J^/'cos^sin^ei(^i-^^+'^) -sin^cos^e5(^^-^i+^)Y (2.12) 
v2 \ 2 2 2 2 J 

In terms of these, the coordinates of the deformed conifold in (12.81) can be written as 



zi = ^{e^^yi + e 2?/i), 2:2 = -^(e2?/2 + e 2^2), 

^3 = ^(e2y3 + e"%3), 2:4 = -^(ei?/4 + e"%4), (2.13) 

where r) is the 'radial coordinate'. The six-dimensional Calabi-Yau metric, dsl, on 
the deformed conifold can be obtained from the Kahler potential given as [41] 

4/3 p 

HV) = ^3]^ rfr/'[sinh(2V) - 2^']'/'. (2.14) 
Here we note that this Kahler potential is derived from the relation 
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= \za\'^ — cosh 77, 



(2.15) 



A=l 



and the metric for the deformed conifold reads as [41] 



dsl 
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(2.16) 



where 



a; = 



u — 



— sin^id(/7i, uj'^ — d9i, 
sin sin ^2C^</'2 + cos ipd92, 



;^ = cos ^l) sin ^2C^</'2 — sin ■0d^2, 

^ — dt/j + cos + COS 92d(p2: 



(2.17) 



and 



X(77) 



(sinh(2?7) — 2?]) 



1/3 



(2.18) 



2^/^ sinh ?7 

The proper radial coordinate which measures the actual distance up the throat in 
the six-dimensional metric is given by 



r{r]) = 



,2/3 



dx 



(2.19) 



VgJo K{x)' 

The warping in this background is induced by the presence of type IIB background 
fluxes including M units of F3 flux through the cycle A and —K units of H-i, flux 
through the cycle B. Due to the presence of three-form fluxes the above background 
emerges as a solution to Einstein's equations and the resulting supergravity back- 
ground is described by a throat with its tip being located at a finite radial coordinate 
riR while at ruv the throat is glued into an unwarped bulk geometry. The warp factor 
is given in terms of the 77 coordinate as [33] 



/iKS = 2(^,Ma')'e-'/-'/(^7). (2.20) 

7(77)= / dx ^^^^^^~^ (sinh x cosh x - xf'^ . (2.21) 
jj^ sinh X 

At small and large radius, the warp factor takes the form: 
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2{gsMa'ye-^/'^ ( 0.5699 - 26-^/3^ J ; rj small 
f^^(lnS + ln|^-lj; ^ large 

In an ISD compactification such as the warped deformed conifold, only the 
complex-structure moduli are stabilized but not the Kahler moduli. Embedding 
the warped throat into the Calabi-Yau sector with Kahler moduli stabilized breaks 
the no-scale structure and induces perturbations around the ISD solution. In the 
noncompact limit (Mpi — )■ oo), these perturbations are sourced by non-linear effects 
including only lASD fluxes. These fluxes satisfy the lASD conditions [17] 

dA = 0, where *6 A = (2.22) 

on general Calabi-Yau cones. The perturbations of $_ around ISD conditions con- 
taining lASD fluxes as the dominant source satisfy the second-order supergravity 
equation of motion [17, 18] 

A(o)<f- = ||AP, (2.23) 
where A(o) is the Laplacian obtained from the unperturbed Calabi-Yau metric g 



(0) 

as before. The solution of Eq. (p.23| ) can always be written as the sum of first order 



homogeneous and second order inhomogeneous solution in the form [17, 18]: 

$_(!/) = $L')(l/) + $L'^(y), (2.24) 
where the homogeneous solution satisfies the Laplace equation 

A^o)^^l\y) = 0. (2.25) 

A particularly simple solution of the Laplace equation on the deformed conifold takes 
the form [32] 

^^}\ri,6) oc (cosh?7sinh?7 — ?7)^''3 cos6'. (2.26) 
The inhomogeneous solution can always be written in the form: 

^^-\y) = 1/ dygiy;y')\A\'iy'), (2.27) 

where 

A^o)Q{y;y')=S{y-y'). (2.28) 

Here and in what follows y denotes a collective internal coordinate consisting of the 
radial and six angular directions, where the later will always be denoted by \E'. On a 
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general Calabi-Yau cone with Kahler form J and holomorphic (3,0) form we may 
turn on (1,2) flux, Ai, and a non-primitive (2,1) flux, A2 [17, 18]: 



Ai = ddhn, 



A2 = dh A J, 



(2.29) 



with /i and /2 being holomorphic functions. The solution (|2.27|) takes the form 
[17, 18]: 

,(2)/ A 9s 



(2.30) 



where /C^*^ is the Kahler metric. For a specific choice of the /j's, Eq. (|2.30|) determines 
the explicit solution of Eq. (|2.23|) up to harmonic terms. 

In the IR region of the throat where 77 is small, the Green's function is that of 
the deformed conifold and takes the form [39] 



Qiy-^y) = Giv,y4; eo) 



32/3 



5Z Yl 1 -^im (1/4,1/4), (2-31) 



28/37^3^8/3 ^ ^ „ 

^=0,^1,..."^=-^ 



where Gq parametrizes the blown up S'^ at = and T^rn, stand for the hypergeo- 
metric functions given by 



2 ' 2 



/45 



= 2^37/1, 



-^0,0 = 1, -^1,1 = 2?/4, 

^1,0 = -v^(l - 42/4^4), -^1-1 = 2y3y4, Ti 
7-ix= -4a/2?/4(1 - 37/4^4), 7-i _i = -4^/2^4(1 - 37/4^4), 



■3 3 

2 ' 2 



4^2/4', 



J-'s 3 

2 ' 2 



4^2^^. 



(2.32) 



In the UV region of the throat where 77 is large, we may introduce another radial 
coordinate r through ~ e^e'' (see Eq. ( [^.15| )) and the Green's function is that of 
the singular conifold given by [13] 



2VAIT4 



1 / r 



X < 



1 / r 



(2.33) 



with the harmonic eigenfunctions 



Yl{^) = Z,„„„«(0i)%,.n2,R(^2)e'"^^'^^+'"^^'^^e5«V', 



(2.34) 
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Zim^A^^i) = Nl{sme,r^ i cot 



2 



R/2 



X 



iFi ( - ii + TTii, + nii, l + TTii-^] sin^ ^ ) , (2.35) 



cot|j X 



2 



„ , R R R r, 0,: . , . 

Fi( -j, + -,l + j, + -,l-m, + -;sm2^ ). (2.36) 



The normalization factors A'^^''^^ relation is given by 

Vti,i / dx[Zj^^rniA^)? / j2,m2,R{y)? = 1- (2-37) 

Jo Jo 
In the above relations 2-^1(0, b, c; d) stands for hypergeometric functions; L is a multi- 
index with the data L = (ji,j2), {^1,^2), R, where ji and j2 are both integers or 
half-integers with mi G {-ji, ■ ■ ■ , ji} and ma G {-ja, ■ ■ ■ , j2}; = 6(ji(ji + 1) + 
^2(j2+l)— -R^/8) denotes the spectrum of the full wave function and the eigenfunctions 
transform under SU{2)i x SU{2)2 as the spin (ji,j2) representation and under the 
U{1)r with charge R; 4 = -2 ± VAl + 4. 

3. Multifield D-brane inflation 

We now embed a mobile D3-brane in the supergravity background described in the 
previous section and analyse its the four- dimensional effective action. In the super- 
gravity background with metric ansatz (|2.2|) the effective action takes the form 



PI 



2 

7dbi = Vl - hgmng''''d^(t)^dr 



/^-'(7dbi-1) + W 



(3.1) 



Here gs is the string coupling, and Mpi is the Planck-mass. The first term in this 
action is the ordinary four-dimensional Einstein-Hilbert action, which arises from 
dimensional reduction of the closed string sector of the ten-dimensional action. The 
second part contains the action that controls the dynamics of the fields, parametrizing 
the position of the brane along the internal coordinates, 0"^. In a strongly warped 
region, 1, the kinetic energy's pre- factor of in Eq. (|3.1|) suppresses it relative 
to V{(f)^) even when the motion is relativistic. 

To study brane inflation, we will take the position of the D3-brane to be ho- 
mogeneous, 0™ = (f)"^{t), and we will consider the four-dimensional metric to be the 
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standard unperturbed FRW metric: 

dsl = g^^dx^'dx" = dt^ - a^{t)dx^, (3.2) 

where a{t) is the scale factor. The variation of the later part of the action ( |3.1| ) with 
respect to the metric produces the energy-momentum tensor, which has the form of 
a perfect fluid with its energy density and pressure given by: 



E = Tsih-'[jj,Bi-l] + V], 
P = T3(/.-M1-7dbi] + ^]- 



(3.3) 
(3.4) 



Here T3 = ((27r)'^a'^) ^ is the D3-brane tension. The full equations of motion are: 



H 



E 



1 d 
dt 



(E + P) 
-T37DBI (7dbi - 1) + 



(3.5) 
(3.6) 



(0) 
In 



kiln 



2h^ 



(3.7) 



where gml. is the unperturbed Calabi-Yau metric, as before. The Planck-mass, Mpi, 
appearing above is bounded by the UV-scale as [32]: 



Ml > ■ 



'pi ^ a J{Vvw), 

where J{r]) = J I{rj) sinh^?7 with I{r]) given by (|2.21|) . 

The D3-brane potential appearing in the above equations of motion is: 

Here the first term is a mass term generated for the canonical infiaton = a 
which by ( ^.10[ ) takes the form: 



(3. 



.2/3 




T3 n dx 

6 Jo 'Kix) 



(3.9) 
(3.10) 



We note that the normalization is strictly dependent on the position of the brane, 
which affects the volume modulus (see e.g. [13-18]), and also on the trajectory of 
the brane, even in the slow roll approximation, due to the infiaton, 0, being a sigma 
model field, [43-45]. The last term in ( p.9| ) arises from perturbations around the 
noncompact (ISD) supergravity solution described by Eq. ( |2.23 ). The homogeneous 
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linearized solution of Eq. (|2.23|) is given by ( p.26| ). The inhomogeneous non-linear 
solution of Eq. (|2.23|) given in its implicit form as ( p. 301 ) equals the nonperturbatively 
generated D3-brane F-term potential in four-dimensional supergravity [18]. The F- 
term potential takes the following form: 



(3.11) 



where {z^} = {p,z^-a = 1,2,3} and Dj^W = dj^W + K,\dj:lC)W with /C and W 
denoting the Kahler potential and superpotential, and = MpJ^ with M^^ bounded 
by (|3.8| ). From ( p.8|) one can see that large r^uv implies large Mpi. In the noncompact 
limit (Mpi — )• oo), the F-term potential ( p.ll|) takes the form 



Vf = JC-^d^Wd^W. (3.12) 
For any holomorphic function W on the conical geometry, we may turn on the flux 

A^gf = d^drWJC^^QQ^p = ddW ■ fi, (3.13) 
where Q is the holomorphic (3, 0)-form. Solving Eq. ( p.23|) for Eq. ( p.l3P gives [18] 



(3.14) 



The Kahler potential, /C, depends on the complex Kahler modulus p = a + ix, and 
on the D3-brane position, {za,Za}, [46] 



/C(;2", p, p) = -3k-' log[p + p - 7A;(z", ^")] = -Sk'' log U{z, p) 



(3.15) 



Here k{z'^, z") is the so-called 'little' Kahler potential of the Calabi-Yau manifold, 
= MpJ^ as before, and 7 is a normalization factor given by 



7 



3M2 



(3.16) 



where (Tq denotes the value of a when the D3-brane is at its stabilized configuration 
[8, 14]. The Kahler metric and its inverse take the form [8] 
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f 1 








U'yk^^ +'y^kak^ J 



Ar 





( U + ^fk^k-^h^ 


k^k^~^ \ 









(3.17) 
(3.18) 



where k^^ = dad^k is the Calabi-Yau metric, and ka = k^^- The superpotential W , 
also depends on the D3-brane positions, {^q,}, and is given by 
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W{p) = Wo + A{z'^)e-''P. 



(3.19) 



The first part of tfie superpotential is tfie Gukov-Vafa perturbative superpotential 
[47]. Tlie second part in tlie superpotential comes from nonperturbative effects 
sourced by moduli stabilizing D7-branes wrapping certain four-cycles in the com- 
pactification and a = 27r/n is a constant with n denoting the number of wrapped 
branes. The infiaton dependence of the superpotential is induced by the interac- 
tion between the inflationary D3-brane and wrapped D7-brane. The displacement of 
the D3-brane in the compactification slightly modifies the supergravity background, 
perturbing the warp factor, h = h^s + ^h, which corrects the warped volume of 
four-cycles from D7-branes. This then corrects the D3-D7 potential inducing the 
backreaction on the mobile D3-brane. The prefactor in the nonperturbative part of 
( p. 19 ) computed from corrections to the warped background takes the form [13]: 



/(O) 



l/n 



(3.20) 



where f{z°') denotes the holomorphic embedding function of D7-branes depending 
on the D3-brane coordinates and Aq is a constant. 

According to this and ( ^.17[ ) - ( p.l^ ), the scalar potential ( |3.11| ) takes the form 
[8, 14] 



K 



2,U{z,pf 

-2 

r 



K 



3Uiz,py 



U{z, p) + -i¥<^k^k-s \W^p\^ - ?>{WW^p + c.c. 



7 



(3.21) 



The first part in Eq. ( p.21|) is the standard KKLT F-term potential, and the second 
part arises exclusively from corrections to the nonperturbative superpotential. By 
considering appropriate formulas for the various terms in the potential (e.g. see [12]), 
it is straightforward to show that the F-term potential ( p. 211) takes the functional 
form Vp = Vf{zi + Zi, \zi\'^,ri,a,x)- To obtain the explicit form of the potential, 
we need to specify an embedding for D7-branes. For the Kuperstein embedding of 
D7-branes [40], we have 



f{zi)=p~zi = 0. (3.22) 

Thus the prefactor of the nonperturbative superpotential, Eq. ( p.20|) , and its deriva- 
tive with respect to the independent coordinates take the form 



A{zi) = Ao[l- 



£i 
p 



l/n 



(3.23) 
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^0 / , Z\ 



A,(zi) = ^ 1-- 5,1. (3.24) 

To obtain a simple trajectory depending only on one harmonic mode, say Q\ = 6, we 
fix the rest of angular directions of T^'^ by imposing the following constrains 



= ± — , = TT, 

2 2' 2 ' 

t^L^£lZ^ = 0, 02 = 0. (3.25) 

Accordingly to this, the coordinates on the deformed conifold read as 



= -e cosh 1 2 j (^2 
Z3 = -iesinh ( | j sin 



Z2 = -ie sinh 1 2 j ( 2 " 



Z4 = +ecosh ( I j sin 



(3.26) 



The imaginary part of the Kahler modulus can be integrated out by 
Accordingly, the four-dimensional supergravity potential is^ 



where 



LP 



9iv,0)\ 



2/n 



\Wo\ e" 



\Ao\g{v,ey/- 



F(r],9) = V/^irSinh^r] 

4 anjjg 



cos - COS - sinh - 



sinhM|) cos^ (^)+^^-^2 



Here we have defined 



(3.27) 



(3.2J 



(3.29) 



U{r],9)=2a{r],9)~jk{r]), 



(3.30) 



""^It is straightforward to check that for 9 = our two-field potential ( p. 28 ) coincides with the 
single field potential derived in [12]. 
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1 H — cosh 



cos 



(3.31) 



Compactifying the warped deformed conifold throat via attaching it to the com- 
pact Calabi-Yau space at some finite radius ruv requires Kahler modulus stabihza- 
tion. This involves integrating out a by the assumption that it evolve adiabatically 
while remaining in its instantaneous minimum cr*(r7, 9) using the following equation: 



9(V^F + Kpiift)(r?,^^) 



da 



0. 



(3.32) 



Here we have added a D-term potential V{rj,a) = D ^puft / U{ri, a Y which is needed 
to uplift KKLT AdS minimum to a dS minimum. The uplifting can be sourced 
by distant anti-D3-branes or wrapped D7-branes. The instantaneous minimum of 
a is denoted by a* including a shift due to its coordinate dependence induced by 
adding a mobile D3-brane to the compactification. The functional form of a^, can be 
determined by the numerical solution of the transcendental equation ( p.32 ), which 
we will solve in the next section. Before turning to numerical computation, we need 
to look at the minimum of the D3-D7 potential which specifies ap and -DupUft and 
constrains the rest of parameters on which these depend. 

The critical value ap of the Kahler modulus before uplifting is determined by 
the condition D„W\r,=o,9=o,ap = 0, or equivalently [2], 



^(l + |a.,),(0.0)V» 



dVp 



•a 



0, 



where g is given by (|3.31|) . We may write this in the form: 



(3.33) 



ap = - log 

a 



l/n- 



From Eq. ( ^281) and Eq. (gjD we note that 

V^F(0,ap) = 



The uplifting parameter is given by 

Kipiift(o, (Tp) 



From this we obtain 



|V^F(0,a^)| 



-Dupiift 



2ap{?, + 2aapy 



with Kplift(0, CXi;') 



-Dupiift 



6s a^K^WQap 
(3 + 2aapY ' 



(3.34) 



(3.35) 



(3.36) 



(3.37) 
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Here 1 < s < 3 to avoid decompactification, a is determined by the choice of n, 
K depends on the UV cut-off, and the value of ap can be derived from Eq. (|3.34| ) 
once the set of parameters {e, /i, n, s, \Aq\, \Wo\ } are specified. We also note that 
uplifting the KKLT AdS minimum to a dS minimum introduces a small shift in the 
stabilized volume, ao = ap + Sa. At the tip, we have [15]: 



aao ~ aap + 



(3.38) 



Here we note that both Eq. (|3.34|) and Eq, 



which give ap and 5a, respectively. 



are derived from the local minimum of the F-term potential ( p. 281) at the tip. The 
critical value am of the Kahler modulus away from the tip and its shift 6am can be 
derived from the global minimum of the F-term potential ( p.28[) . By computing the 
first and the second derivatives of the F-term potential Vp with respect to cr, it is 
straightforward to show that the global minimum of Vp requires: 



ar. 



l/n 1^0 1 



aU aU + 2 

1 + — + 2aF——— 

3 all + A 



(3.39) 



Under the assuption that Vd only shifts the minimum by a small amount gives: 



V'{am + 6a) = Vi,{am + 6a) + V^{am + 6a) ~ Vi,{am) + 6aV;^{am) = 0(3.40) 

Vl^{am) UmVp{am)V"{am)' ^ ■ ^ 

We would like to remark here that in most of the brane inflation models in the 
literature (e.g. see [12, 15]) the perturbative expansion of the potential is analysed 
in either the UV or the IR region. In these models e and fi are considered to have 
similar order magnitude. However, we note that a consistent expansion along the 
entire throat including both the UV and IR regions requires a large hierarchy between 
e and /i. In particular, the final piece of our potential term F{ri, 9) given by Eq. (|3.29|) 
scales as e~^/^(e//i)^. Unless the hierarchy between e and fi is large this term will 
dominate the potential and destabilize the vacuum expansion. Moreover, the key 
point here is that in our DBI brane inflation set up the full potential contains a large 
leading order mass term for the radial coordinate and in order to keep the adiabatic 
approximation ( p.32|) valid the hierarchy between e and fi has to be large, so that 



the mass generated for the Kahler modulus is much larger than the Hubble rate (see 
Section 4). 

As mentioned above, apart from the shift induced by uplifting, 6a, the addition 
of a mobile D3-brane in the compactification induces a further shift in the Kahler 
modulus that depends on the coordinates of the brane, a^{rj,6). Thus the nonper- 
turbatively generated D3-D7 potential about a dS minimum takes the form: 
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vp + vo= i °' 1 .... 



uplift 



Hence the total potential takes the form: 



(3.42) 



T3 yjj^l Vilf + C2K{r]) sinh ?7 cos -i?] + Vq 
+Vf + Vd. 



(3.43) 



Here I^upUft is given by (p.37|) , the constant Vq is chosen so that the global minimum 
of is = 0, and C2 is an arbitrary constant. Because C2 multiphes a solution to a 
free Laplace equation, it is not fixed. For a self-consistent expansion, we expect C2 
to be smaller, or of a magnitude comparable with other terms in the potential. 

Finally, to derive the explicit form of the D3-brane equations of motion on the 
deformed conifold for the simplest case including only one angular direction for the 
potential ( p.43 ) first note that a simple round metric on the deformed conifold 



can be obtained from (2.16) as: 



ds^ = A{r])d7]^ + B{7])de\ (3.44) 

where 

A{ri) = B{^) = ^M[cosh(r//2) + sinh(r//2)]. (3.45) 



The brane equations of motion can then be derived upon cross elimination from (|3r 
in the following form: 



^^-^ + -^^'(1 - 7DBi) + 7T7W(7dbi - 1)' 



7DB1T3 jBBiAT'^ 



3 
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me 

7dbi 



;i - 7dbi)^^7 



h' 



Or]— + Wr] ^ 

B Tdbi^s 



-fl - hB9^] 



-fl 



7dbi) 



1 



7dbi'^ 



2hB 



h 



(3.47) 



Here we note that in the absence of non-hnear corrections including the contribution 
of the D3-D7 potential and perturbations of the warp factor, the potential ( p.43| ) 
and Eqs. ( |3.46|) - ( p.47|) reduce to the D3-brane potential and equations of motion in- 
cluding only linearized corrections studied in [32] . We also note that in the slow-roll 
regime 7dbi — 1 and the perturbations of the warp factor have no effect. However, 
since we are neither slow-rolling nor restricting ourselves to the linearized case, we 
need to take into account the perturbations of the warp factor which amount super- 
potential corrections in the D3-D7 potential. For the warped deformed conifold, the 
perturbations of the warp factor in the tip region can be expressed as [39] 



5h 



16n ■ 3^/^g,p{a'f 



(2e)8/3 ■ T] 



1 + 2V2yi + 6yl + %V2y, 



OA . fe2\ 

COS I — I sm I — I e 



sm I — I cos I — I e 



with (3.48) 



. (3.49) 



Here Qiji^yi) is the Green's function (expanded in the eigenfunctions of the Lapla- 
cian) on the deformed conifold given by (|2.31|) , ?/4 comes from ( p.l2[) and p specifies 
the number of mobile D3-branes, which we take to be p = 1. By the first angular 



condition in (p.25|) along an S' 



0) for relations (|3.48|) and ( p.49| ) we have 



5h 



167r ■ 32/3^,p(a') 
(2e)8/3 ■ r] 



1 . fe 

,4 = -^sm 



1 + 2722/4 + Qyl + Sv^yl - ■ ■ 



with 



(3.50) 
(3.51) 



In the off-tip region, where r is large, we may use the very first line in ( |3.48D together 
with relations (|2.33|) - (|2.37|) for the quantum numbers ji =32 = R/2 = 1/2 and 
mi = m2 = 1/2, and obtain 5h in the form 



,/2 



5h 



lOgspa 4 

n cos - 

7r2 2 



9 ( 1 



with 



(3.52) 
(3.53) 
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Figure 1: The behaviour of ap for choice of compactification parameters 
(4.1), and the behaviour of for choice of parameters (4.1) with 6 = 

(sohd), 7r/2 (tiny-dashed) and vr (large-dashed). 




Figure 2: The behaviour of the function F and 5a for the choice of parameters (4.1) and 




t igure 3: The Hubble rate, if^ ~ V/3M^^ , and the a mass squared, for the choice 
of parameters (4.1) and 6 = Tr. 



4. Inflationary solutions 

In order to integrate the full brane equations of motion, we first need to specify a suit- 
able choice of parameters and then compute the real part of the Kahler modulus that 
appears in th e brane equations of motion. To choose a reasonable set of compactifi- 
cation parameters, we note the following points. Firstly, we require a large hierarchy 
between Aq and Wq to guarantee large ap which ensures suppressed a'-corrections. 
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Secondly, we also need a large hierarchy between e and /i in addition to the Aq/Wq 
hierarchy to guarantee a valid perturbative expansion, 6a ^ 0. When both these 
hierarchies are turned on, a^, can be computed within the adiabatic approximation 
from Eq. (|3.32|) . We also remark from the literature that choosing a large value of 
the UV-scale, rjuy, sets a large value for the Planck-mass (e.g. see [32]) in which 
case curvature corrections may be omitted and non-linear corrections are dominated 
by lASD fluxes sourced by moduli stabilizing wrapped D7-branes whose number is 
given by n > 1. Furthermore, we remark that the supergravity solution requires 
large ggM and the value of s has to be chosen within the range 1 < s < 3 to ensure 
a small positive cosmological constant and to avoid runaway decompactification. 

In line with the above requirements we choose in our numerical analysis the 
following specific set of compactification parameters: 



T]uY = 15, n = 2, s = 2, e = 0.001 

gsM = 100, /i = 5, Wo = 29, = 25 x 10^^ (4.1) 



Inspection of Eq. (|3.34|) and Eq. (|3.39|) shows that for the choice of parameters (|4.1j ) 



ap and am scale to large values (see Fig.|^); inspection of Eq. ( p.29| ) and Eq. (|3.41 



shows that for the choice of parameters ([4.1|) and 9 = it the functions 6a and F scale 
to reasonably small values on the entire warped deformed conifold (see Fig.|^); for 
the choice of parameters ( |4.1j ) the mass squared of the Kahler modulus, = V„a, 
is much larger than the approximate squared Hubble rate, ~ V/SM^^ (see Fig.|^), 
which guarantees an adiabatic expansion [15, 48]. We note here that by decreasing 
these hierarchies much below their values ([l.l]) makes and scale similarly 



which invalidates the adiabatic approximation, in addition 6a and F become large 
which invalidates the perturbative expansion. 

In the literature of brane inflation, the most common way of computing the 
Kahler modulus, a, is to adopt a semi - analytic approach with the assumption that 
a evolves adiabatically (e.g. see [14]). In this approach, a in U{ri, a) is set to its large 
fixed value a^, and (|3.32|) is treated as an equation in the variable exp[— acr*(?7)]. The 
value of a^: is then obtained from Eq. ( p.32| ) by expanding in specific conical regions 
including either the IR or the UV regions where the canonical inflaton (|3.10|) is 
given by its small or large rj limit, respectively. This may be a good approximation 
but it gives only a qualitative understanding. Here we compute a^, by solving the 
transcendental equation ( |3.32| ) numerically to obtain the exact value of a^, on the 
entire supergravity background. For the choice of parameters ( [4.1[ ) we show a^: in 
Fig.|. 

For the choice of parameters (f4.1|) and the numerically computed Kahler mod- 
ulus, we integrated the full D3-brane equations of motion, Eqs. ( p.5[ ) - ( p.7| ) with 
Eq. ( prrp given by Eqs. (|3.46| ) - (|3.47| ), and our inflationary solution is displayed in 
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Fig.^. The solution describes spiral brane motion at high speed in the warped 
throat region of the compact Calabi-Yau space containing holomorphically embed- 
ded wrapped D7-branes involved in (Kahler) moduli stabilization. The conserved 
angular momentum is lifted by harmonic dependent corrections from linearized as 
well as non-linear perturbations including contributions from the D3-D7 potential 
and corrections of the warp factor. The brane accelerates along the radial and angu- 
lar directions as it falls down the throat from the UV end where it is attached to the 
compact Calabi-Yau space. Inflation ends when the brane reaches the IR location 
where the throat smoothly closes off. For the choice of parameters and the numer- 
ically computed Kahler modulus, we integrated the brane equations of motion and 
found that the inflationary solution is quite robust against harmonic dependent cor- 
rections from the D3-D7 potential and corrections of the warp factor. In particular, 
we found (as displayed in Fig. |) that harmonic dependent corrections induced by 
the D3-D7 potential and perturbations of the warp factor do not have the effect of 
increasing the number of e-foldings and decreasing the 7DBi-factor. This result differs 
from our previous results [32] in which harmonic dependent correction to brane mo- 
tion from linearized perturbations of the supergravity solution increased the number 
of e-foldings compared to the number of e-foldings produced by brane motion with 
conserved angular momentum (spinfiation) where no supergravity corrections and 
hence no harmonic dependence in brane motion is present. 

We repeated the above computation for a various choices of initial conditions 
and compactification parameters and our findings are summarized as follows. 

• Decreasing e below its considered value while keeping other parameters fixed tends 
to flatten the functional form of cr* and changes its overall scale insignificantly. Also, 
the decrease in e leaves ap unchanged. 

• Decreasing the value of /i while keeping other parameters fixed slightly decreases 
the value of ap and changes the scale of a^, by a minimal amount. Note that here we 
decrease n by an amount, so that the hierarchy between e and fi still remains large. 

• Increasing/decreasing n while keeping other parameters fixed strongly impacts the 
ap and the scale of cr*. Also, increasing n induces fluctuations in a* at large rj. 

• Increasing/decreasing the hierarchy between Aq and Wq while keeping other pa- 
rameters fixed increases/decreases the value of ap and slightly changes the scale of a^, 
but leaves its overall shape unchanged. Here again we do not decrease the hierarchy 
between Aq and Wq too much. Also, increasing the hierarchy between Aq and Wq by 
taking Wo much smaller than its considered value increases ap and suppresses the 
uplifting contribution. This slightly changes the form of a^, near the origin. 

• Changing the set of compactification parameters and initial conditions indicates 
that the number of e-foldings produced by spinfiation including non-linear harmonic 
dependent corrections depends more on the subset of compactification parameters 
{e, gsM,7]ijy} than on the initial conditions and is insensitive to the choice of the 
remaning subset of compactification parameters {n, s, fi, Aq,Wq}. Comparing with 
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Figure 4: The Kahler modulus, a^:{ri,9), on the entire supergravity background with the 
choice of compactification parameters (4.1). 



the number of e-foldings generated by spinflation including only linearized harmonic 
dependent corrections each time when varying the set of parameters and initial con- 
ditions shows no differece in the number of e-foldings (as Fig. ^). 
• The 7DBi-factor produced by spinflation including non-linear harmonic dependent 
corrections is insensitive to the choice of compactification parameters and initial con- 
ditions. Comparing with the 7DBi-factor produced by spinflation with only linearized 
harmonic dependent corrections each time when varying the set of parameters and 
initial conditions shows no differece in the 7DBi-factors (as Fig. H). 

The above findings show that the inflationary solution is quite robust against 
harmonic dependent corrections from the D3-D7 potential and perturbations of the 
warp factor not just for the specific choice of compactification parameters ([4.1|) but 
for a very large set of consistent parameters. 



5. Summary and conclusions 

In this paper we studied brane inflation in a warped string compactification incor- 
porating the effects of moduli stabilization and backreaction from UV-deformations 
of the warped throat geometry. The focus of our paper was on DBI brane inflation 
in the warped deformed conifold with a UV/IR consistent perturbative expansion 
around the noncompact ISD solution. The perturbations were dominated by lASD 
fluxes sourced by moduli stabilizing wrapped D7-branes. 
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t 



Figure 5: The number of e-foldings, A^, and the 7-factor, 7dbI) with (gray-dashed) and 
without (black-sohd) non-hnear harmonic dependent corrections for the choice of compact- 
ification parameters (4.1). 



We computed the D3-brane potential on the entire deformed conifold includ- 
ing non-linear corrections from the flux induced potential in ten-dimensional su- 
pergravity which equals the nonperturbatively generated D3-D7 brane potential in 
four-dimensional supergravity. For a simple choice of a trajectory on the deformed 
conifold, we integrated out the Kahler modulus and reduced the D3-brane potential 
to a simple two-field potential depending on radial and harmonic directions of the 
deformed conifold. We integrated out the Kahler modulus by full numerical compu- 
tation determining its exact functional form on the entire supergravity background 
including both the IR and UV regions. We found that a UV/IR consistent pcrtur- 
bative expansion in the supergravity potential with the Kahler modulus integrated 
out within the adiabatic approach in DBI inflation requires certain hierarchies of 
scales that determine the set of compactification parameters different from those in 
slow-roll models. 
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For the consistent choice of parameters and the numerically computed Kahler 
modulus, we integrated the D3-brane equations of motion in the warped deformed 
conifold with harmonic dependence from the D3-brane potential and perturbations of 
the warp factor. We found that our numerical solutions are quite robust against non- 
linear perturbations including harmonic dependent corrections from perturbations of 
the warp factor and the D3-D7 brane potential. In particular, we found that harmonic 
dependent corrections from the D3-D7 potential and perturbations of the warp factor 
do not have the effect of increasing the number of e-foldings and decreasing the 
7DBi-factor. We therefore conclude that the most leading order harmonic dependent 
correction to brane (spin) inflation comes from the linearized corrections analysed in 
[32] with the level of non-Gaussianity remaining large. 

Our analysis can be extended in several ways. One direction would be to con- 
sider different embedding functions for D7-branes with different trajectories on the 
deformed conifold and see how this may affect the inflationary solutions. Despite 
the fact that different embedding functions for D7-branes and different trajectories 
on the deformed conifold modify the functional form of the supergravity potential, 
the Kahler modulus and perturbations of the warp factor, we expect this to have a 
subdominant effect on the number of e-foldings and the 7DBi-factor. It would also 
be interesting to consider the possibility of a dynamical Kahler modulus (instead 
of stabilized) [48], and integrate the brane equations of motion for a less restricted 
choice of parameters. Since taking the Kahler modulus field dynamical adds to the 
number of brane equations of motion, we expect this to have a less trivial impact on 
the inflationary solutions, though we do not expect this to decrease the 7DBi-factor 
by an appreciable amount. 

The other, perhaps more interesting way of extending our analysis is to consider 
further corrections to the infiaton action and analyse in detail the effects of cosmo- 
logical perturbation theory. One particular correction comes from the contribution of 
the flux induced potential of harmonic type [18]. The flux induced D3-D7 potential 
that we considered in our inflationary analysis came from the holomorphic solution 
of the noncompact supergravity equation of motion which described non-linear per- 
turbations around the ISD solution. Since the general solution should be harmonic 
rather than just holomorphic it would be necessary to include such harmonic contri- 
butions which have not been computed for the deformed conifold to date. 

Another correction to the inflaton action arises from departures of the noncom- 
pact limit [18] considered in this paper. Despite the fact that our Planck-mass was 
set large by the UV-scale we considered, it would be interesting to consider possible 
departures from the noncompact limit and compute further contributions to the D3- 
brane potential from coupling to curvature corrections which also induce harmonic 
dependence in brane motion [18]. In particular, coupling to the Ricci-scalar intro- 
duces a non-minimal coupling in the DBI action which corrects the 7DBi-factor and 
may have the capacity of decreasing the level of non-Gaussianity significantly [49]. 
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It would be interesting to confirm this result in the concrete supergravity set up con- 
sidered in this paper and investigate its implications for cosmological perturbations 
in more detail [50] in our framework. 

Finally, it would be important to confirm whether in our supergravity set up mul- 
tifield effects (e.g. from phase transition) induced by an instability along harmonic 
directions [51] do have the capacity to evade stringent constraints in cosmological 
perturbations for single field infiationary models. Moreover, it would be very in- 
teresting, though formidable, to consider our DBI brane inflation model along a 
trajectory on the deformed conifold depending on all six directions and extract the 
full multifield effects for cosmological perturbations in a UV/IR consistent expansion 
and make contact with some of the results obtained in [52] by taking the singular 
conifold limit. We shall leave the investigation of these for the future. 
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